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Abstract: We explore general SUSY GUT models with exact third-generation Yukawa 
unification, but where the requirement of universal soft terms at the GUT scale is relaxed. 
We consider the scenario in which the breaking of universality inherits from the Yukawa 
couplings, i.e. is of minimal flavor violating (MFV) type. In particular, the MFV principle 
allows for a splitting between the up- type and the down-type soft trilinear couplings. 
We explore the viability of this trilinear splitting scenario by means of a fitting procedure to 
electroweak observables, quark masses as well as flavor-changing neutral current processes. 
Phenomenological viability singles out one main scenario. This scenario is characterized 
by a sizable splitting between the trilinear soft terms and a large \x term. Remarkably, 
this scenario does not invoke a partial decoupling of the sparticle spectrum, as in the case 
of universal soft terms, but instead it requires part of the spectrum, notably the lightest 
stop, the gluino and the lightest charginos and neutralinos to be very close to the current 
experimental limits. The above mechanism is mostly triggered by a non-trivial interplay 
between the requirements of negative, sizable SUSY threshold corrections to m& and an 
instead negligible modification of the B — > X s ^f decay rate, in presence of various other 
constraints, most notably a successful EWSB and a not too large BR(f? s — ► 
We present a model-building interpretation of our discussed scenario and emphasize the 
crucial role of SUSY spectrum determinations at the LHC for either falsifying Yukawa 
unification or else providing important hints on the mechanism of SUSY breaking at work. 
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1. Introduction 

It remains a striking fact that the Standard Model (SM) gauge couplings, measured at low 
energies, evolve to a single value at a high scale Mq, provided that, above the electroweak 
scale, the SM is assumed to become supersymmetric. This fact may be an accident or 
may not be. In support of the second possibility are a few further remarkable features of 
this observed gauge coupling unification: it is very weakly dependent on the details of the 
supersymmetric spectrum, hence (presumably) a robust consequence of the assumed gauge 
symmetry and of supersymmetry (SUSY); furthermore, Mq ~ 3 x 10 16 GeV happens to be 
at just the right place, namely above the region where proton decay at an unacceptable rate 
is generic, and below Mpi anc k, where unavoidably large gravitational effects would make 
the calculation unrealistic. These considerations, together with the possibility to address 
various structural questions unanswered within the SM, warrant the decade-long interest 
in grand unified theories (GUTs). However, concerning further tests beyond that of gauge 
coupling unification, the possibility of general conclusions has been hampered by the larger 
number of model assumptions needed in each case. A prototype example is the already 
mentioned proton decay [1, 2, 3, 4, 5, 6, 7, 8]. 

An alternative way to test SUSY GUTs is by exploring the consequences of the generic 
expectation of low-energy SUSY. In general, the question of the predicted pattern of SUSY 
masses and mixings maps onto the question of the mechanism of SUSY breaking and of the 
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form Yukawa couplings assume at the high scale. In turn, the latter issues usually require 
strong theoretical assumptions, that however may be justified, besides their possible appeal, 
if they lead to sufficiently sharp predictions. 

A very elegant assumption, potentially testable in the SM fermion masses and mixings, 
is that of Yukawa unification (YU) at the GUT scale [9, 10, 11, 12, 13, 14, 15, 16]. It is 
motivated by the fact that, due to the higher degree of symmetry, matter fields must sit in 
appropriate representations of the gauge group, thereby sharing a common Yukawa cou- 
pling. Since this simple picture can be spoiled by e.g. the presence of higher-dimensional 
interactions, the crucial question is whether YU may leave any low-energy remnant at all. 
While for the light fermion generations this is definitely not the case, 1 for the third genera- 
tion it remains an open and appealing possibility, selecting the group SO(10) as the poten- 
tially most predictive case. In this case, verification of the YU hypothesis would amount 
to the important conclusion that all the dimension-4 interactions involving 3rd-generation 
fermions and/or scalars originate from just the I63 I63 10h structure. More shaky are 
symmetry assumptions, e.g. universality, on the soft-breaking terms, which correspond to 
operators of dimension less than 4, and are hence unlikely to preserve information about 
the symmetries inherent to the UV theory completion. The tacit motivation here is just 
one of computational simplicity. 

In ref. [19] the viability of the hypothesis olt — b — r Yukawa unification in SUSY GUTs 
was studied, assuming that soft-breaking terms for sfermions and gauginos are universal at 
the GUT scale. It was found that this hypothesis is challenged by the constraints imposed 
on the parameter space by FCNC processes, unless decoupling of the squark spectrum is 
invoked, thereby pushing the lightest squark well above 1 TeV. Under the same universality 
hypotheses, a viable alternative to decoupling has been found to be a moderate breaking 
of t — b unification while keeping b — r unification, or equivalently a parametric departure 
of tan j3 from the value implied by exact t — b — r YU. 

The conclusions of ref. [19] hold, we repeat, under the assumption, very common in the 
literature, of GUT-scale universalities for soft terms in the sfermion as well as in the gaugino 
sector. 2 An interesting question is then whether departures from GUT-scale universalities 
exist, that on the one hand allow Yukawa-unified SUSY GUTs to successfully withstand 
the combined constraints mentioned above without decoupling of the SUSY spectrum, and 
that on the other hand can be substantiated with a plausible SUSY-breaking mechanism. 

In this respect, two simple and well-motivated scenarios of GUT-scale non-universalities 
emerge [20]: 

• Non- universal gaugino masses (NUGM), 

• Non-universal scalar soft terms of minimally flavor-violating (MFV) form, i.e. inher- 
iting from the SM Yukawa couplings. 

1 Yukawa ratios different from 1, which are predicted in some GUT models, can however be phenomeno- 
logically viable [17, 18]. 

2 On the other hand, soft terms for the Higgs scalars are allowed, and actually required, to be split from 
each other. 
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These possibilities are not exclusive to each other - e.g. there could even be a single i^-term 
SUSY breaking causing non-universalities of both the above mentioned types at the same 
time [20]. However, to simplify matters, we will keep these two possibilities separate in 
the rest of the present discussion. In particular, since a detailed analysis of the NUGM 
scenario can be found in Baer et at, ref. [21], and a dedicated study for the case of YU in 
[22], we will focus on the scenario of MFV soft terms. 

Concretely, the purpose of the present paper is to address the following questions 

1. whether relaxing GUT-scale universalities in favor of MFV soft terms allows, in the 
context of SUSY GUTs with YU, to recover phenomenological viability without in- 
voking decoupling of the sfermion sector; 

2. whether the experimental constraints used to address point 1 provide enough infor- 
mation to single out specific regions in the general parameterization, to be discussed 
below, that soft terms respect when they are of the MFV form. 

The question of phenomenological viability, point 1, will be addressed by contrasting 
our class of models with established data on EW observables and flavor-changing neutral 
current (FCNC) decays. Our conclusions will be assessed through a fitting procedure, simi- 
lar to that of refs. [23, 19] , which has the advantage of being manifestly reparameterization- 
invariant. The details of this procedure will be presented in section 3. 

Turning to point 2, the corresponding question is relevant in connection with the search 
of a plausible mechanism of SUSY breaking able to substantiate the pattern of soft terms 
emerging from point 1. As anticipated at the end of the next section and detailed in our 
numerical analysis, section 4, low-energy constraints are indeed powerful enough for a clear 
pattern of soft terms to emerge. A concrete example of a SUSY-breaking scenario where 
this pattern naturally emerges will be discussed in section 5. 

Our results provide a concrete example where, under what we consider very plausible 
(and to our knowledge previously unexplored) assumptions for soft terms, enough remnant 
information on the high-energy symmetries survives at low energies for these symmetries 
to be reconstructible. Our results also illustrate the crucial role that measurements of the 
lightest part of the SUSY spectrum play in this reconstruction program. 

2. Yukawa unification and the MFV principle 

In order to introduce our problem of interest, let us first review briefly the case of Yukawa- 
unified SUSY GUTs where soft terms at the GUT scale are parameterized in terms of a 
universal soft mass mi6, a universal trilinear coupling Aq and a universal gaugino mass 

The assumption of unification of the third generation Yukawa couplings at the scale 
Mq ~ 3 x 10 16 GeV gives rise to two relations between the top quark, bottom quark and tau 
lepton masses. These relations depend on the renormalization group (RG) evolution, which 
is governed by gauge and Yukawa couplings only, and on weak scale threshold corrections, 
which depend on the soft SUSY breaking parameters, or, equivalently, on the sparticle 
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spectrum and mixings. These threshold corrections are most important for the b quark 
due to non-holomorphic contributions enhanced by tan/3 [24]. Neglecting these threshold 
corrections, requiring the tau lepton mass to equal its observed value of 1.777 GeV and 
choosing tan (3 ~ 50 in order to reproduce the measured top quark mass [25], the running b 
quark mass ra^m;,) would be predicted to be 4.5 GeV, as opposed to the precisely measured 
experimental figure of (4.20 ± 0.07) GeV [26]. This illustrates the necessity of considering 
regions of SUSY parameter space where the overall threshold corrections to m\, are negative 
[27, 28]. We will now discuss for which choices of GUT-scale parameters this is the case. 

There are two dominant contributions to the threshold corrections to m&, one arising 
from gluino-sbottom loops and one from chargino-stop loops. Considering only these two 
dominant contributions, 3 the running b quark mass at the decoupling scale can be written 
in terms of the running b quark Yukawa coupling and the threshold corrections as [24] 



= ^|cos/3(l + A3 + A^) , (2.1) 



where 



A g = ^^Aitan/3m 5 /(m Si ,mg 2 ,m § ) , (2.2) 

and the loop function, which (for positive arguments) is strictly positive and has dimensions 
of inverse mass-squared, is 

xylnx/y + yzlny/z + zxlnz/x 

I(a, b, c) = - r . (2.4) 

(x - y)(y - z){z - x) 

Taking the /i parameter to be positive, as is indicated by the muon {g — 2) anomaly, 
Ag leads to a positive correction of nib, while the sign of is given by the sign of the 
stop trilinear parameter At. The assumption of fi > will be kept throughout this paper 
as well. To fulfill the condition A g +Ay, < 0, it is therefore necessary to have At large and 
negative. 

Assuming a universal trilinear coupling Aq and a universal gaugino mass mjrt at the 
GUT scale, the low-energy value of A t is given by 

A t » -2.0 mi/a + 0.23 A , (2.5) 

hence large negative At requires large negative Aq or large m\i<i. However, the latter 
possibility is precluded by the fact that the running gluino mass is given by 

m~ g = M 3 ^ 2.6 m 1/2 , (2.6) 

3 In the numerical analysis of section 4, we take into account all contributions to the threshold corrections. 
4 Our sign convention for At is such that the off-diagonal entry of the tree-level stop mass matrix reads 
mt{A t — MCot/3). 
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hence small mi /2 is required to suppress A^. 5 

Even with \At\ 3> Tn\/2, the gluino contribution to mj is still competitive with the 
chargino contribution in a large portion of MSSM parameter space. This is why an ad- 
ditional suppression is necessary, which can be achieved by a large stop-sbottom mass 
splitting, so regions of parameter space where the stop is the lightest sfermion are pre- 
ferred. 

All the features described above can be realized in the framework of the NUHM, the 
MSSM with non-universal Higgs mass parameters, and indeed in refs. [27, 28] (see also 
[29]) the region with 

-A w 2 mi 6 , H, mi/2 < mi6, (2.7) 

was found to allow successful YU. Within the parameter space of eq. (2.7), successful 
EWSB requires Higgs mass-squared parameters with the pattern m\ G < TTVjj^ < vr? H . 

However, relations (2.7), together with the large value of tan/3 ~ 50 required for YU, 
have an important impact on the SUSY spectrum and on the predictions for FCNCs, in 
particular on those B- physics decay modes that are especially sensitive to large tan/3 and 
to the large trilinear coupling At implied by relations (2.7). Specifically, the decay modes 
that turn out to have the strongest impact are B s — ► B — > X s j and B — ► X s £ + £~ . 

For example, in the parameter space (2.7), one typically has to face a substantial enhance- 
ment of BR(j5 s — ► /i + /i~) and huge destructive interference from chargino contributions in 
BR{B X a i). 

In ref. [19] the non-trivial interplay among these observables and the bottom mass has 
been studied extensively through a fitting procedure. The main conclusions were: 

(a) that a generic SUSY GUT with exact YU and GUT-scale universalities for sfermions 
and gauginos is phenomenologically viable only by advocating partial decoupling of 
the sfermion sector, the lightest mass exceeding 1 TeV; 

(b) that phenomenological viability can be recovered without decoupling by relaxing 
t — b — t unification to b — r unification, equivalent to a parametric departure of 
tan/3 from the value implied by exact YU. This solution is non-trivial since, while 
the FCNC constraints prefer lower values of tan /3, a successful prediction of in 
YU requires instead a value of tan/3 very close to 50 [30, 31]. Indeed, a compromise 
solution between the two classes of constraints has been found to exist only for the 
narrow range 46 < tan/3 < 48, implying that the breaking of t — b YU must be 
moderate, in the range 10 - 20%. 

As stated in the introduction, our aim here is to address instead whether the problem 
in point (a) can be reconciled with the assumption of exact t — b — r Yukawa unification 
by relaxing instead the strong (and theoretically poorly justified) hypothesis of GUT-scale 
universalities in soft terms. 

5 This tension can be relieved by allowing non-universal gaugino masses at the GUT scale. In that case 
eq. (2.5) generalizes to A t « -0.2M 2 - I.8M3 + 0.2SA . As mentioned in the introduction, this NUGM 
scenario can lead to viable models of YU, but will not be considered in the following. 
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We will focus in this paper on non-universal soft terms for scalars satisfying the prin- 
ciple of minimal flavor violation (MFV) [32, 33], corresponding to the assumption that 
the threshold at which the flavor symmetry is broken lies above the scale at which the 
soft terms are specified, and that the only spurions of the broken flavor symmetry are 
the Yukawa couplings of the corresponding SM interactions, so that the soft terms flavor 
structure must inherit from the SM Yukawas themselves. Then, restricting to the partners 
of the quark sector, soft scalar mass and trilinear terms have the form [20, 32] 

m 2 Q (l + c u Q YuYi + c d Q Y D Y D + 0{Y^ D )) , 

mU^ + c^Yu + OiYS)) , 

m 2 D {l + c d D YlY D + 0{Y D )) , 
A v = A U Y U (1 + 0(Y 2 )) , 

A D = A D Y D {1 + 0{Y 2 )) , (2.8) 

where the fn and A parameters have mass dimension 1 and the c parameters are real, O(l) 
numbers. 

The hypothesis of exact YU - and use of the hierarchical structure of Yukawa cou- 
plings - allows to drastically simplify expansions (2.8). The SUSY-breaking terms in eq. 
(2.8) assume in fact the approximate pattern 



m Q 



m Q,U,D 



( m Q,U,D \ 

m Q,U,D 
m 2 Q 

,U,D + Am Q,U,D/ 



(2.9) 



Au{D) 



/ 

(2.10) 

\ y t (b)A u{D ) 



i.e. they can be taken as diagonal and split between the third and the first two generations. 6 
This approximation is phenomenologically valid up to terms of the order (Y^ D )ij/y^ 3 with 
i,j 3 and 2/33 the common value of the GUT-scale Yukawa coupling for the third gen- 
eration fermions. In particular, this approximation is sufficient to reproduce the most 
important features of the low-energy sparticle spectrum, while the neglected off-diagonal 
terms are of subleading importance for the spectrum, and also of subleading importance 
for the success of YU. 

Therefore, in the instance where the SM flavor symmetry group is broken minimally 
at a scale higher than Mgut ; the hierarchy of Yukawa couplings and the assumption of 
YU allow to parameterize the GUT-scale soft SUSY-breaking terms in the squark sector 
in a generic way with 6 real parameters for bilinear soft terms and 2 complex parameters 
for trilinear soft terms. 

Given this general parameterization of squark soft terms, three scenarios of non- 
universalities compatible with the MFV principle suggest themselves: 



°To establish contact with A t: b, used before in the text, we note that (^4^0)33 = yt,bA t 
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1. generational bilinear splitting, Amg v D ^ , 

2. up-down bilinear splitting, Wlq ^ mfj ^ m 2 D , 

3. up-down trilinear splitting, Ajj ^ Ajj . 

Concerning the first of these possibilities, it is interesting to note that the viability of 
YU, as described at the beginning of this section, is essentially determined by quantities 
related to the third squark generation, such as the stop and sbottom spectrum; further- 
more, conditions (2.7), which were found to be favourable to YU, naturally lead to a large 
hierarchy between light third generation squarks and heavy first and second generation 
squarks, i.e. to an inverse scalar mass hierarchy [34]. Thus, changing the masses of first 
and second generation squarks by means of generational bilinear splitting is not expected 
to have a strong impact on the success of YU. Indeed, our initial numerical explorations 
of YU in this scenario pointed to preferred values of Arriq v D ~ 0. Hence, we will not 
consider any generational splitting in the following. We note however that it might help 
accommodate additional constraints, like the dark matter relic density, while not upsetting 
the mechanism ensuring the success of YU [35] . 

Scenario 2. instead can have a more profound impact on YU, since the different 
GUT-scale values for the up-type and down-type squarks can lead to a larger hierarchy 
between the stop and sbottom masses at low energies than is possible in the universal case. 
This would in turn allow a suppression of the unwanted positive gluino contributions to 771&. 
Indeed, such scenario has been studied before in the context of b — r unification [36, 37, 38]. 
Whether it is possible to accommodate full t — b — T YU and to satisfy all FCNC constraints 
within this scenario would require a dedicated analysis, which we however leave to a future 
study. 

For the remainder of this work, we will thus concentrate on the study of the trilinear 
splitting scenario, where by definition a splitting between the up-type and down-type tri- 
linear couplings is assumed, whereas sfermion bilinears are still taken as universal. That 
is, our assumptions for the soft terms of squarks and sleptons at the GUT scale are 

m Q,u,D,L,E = mle 1 » ( 2 - n ) 
A U = A U Y U , A D = A D Y D , A L = A D Y L . (2.12) 

3. Procedure 

Our problem of interest is to study the MFV scenario with soft terms defined by eqs. 
(2.11)-(2.12) and a positive fi parameter in the framework of a generic SUSY GUT with 
exact t — b—T YU. The present section is devoted to the discussion of the procedure adopted 
in our numerical analysis. 

The parameter space of our considered class of models is constrained through a fitting 
procedure against low-energy observables, that are reported in tables 1 and 2 along with 
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Observable 


Value (cr exp ) 


Observable 


Lower Bound 


M w 


80.398(25) 




114.4 


M z 


91.1876(21) 


M x+ 


104 


1 r&ri 
1U G M 


l.lDDO ( {I) 


M t 


wo. i 


l/a cm 


137.036(0) 






a s (M z ) 


0.1176(20) 






M t 


173.1(1.3) 






m b (m b ) 


4.20(7) 






M T 


1.777(0) 







Table 1: Flavor conserving observables [26, 25] used in the fit. Dimensionful quantities are ex- 
pressed in powers of GeV. 



Observable 




Value((Texp)(o-theo) 


Ref. 


AM s /AM d 




35.1(0.4)(3.6) 


[39, 40] 


10 4 BR(B - 




3.52(25)(46) 


[39] 


10 6 BR(B - 


> X s £+£~) 


1.60(51)(40) 


[41, 42] 


10 4 BR(B + 




1.40(40)(26) 


[26] 


BR(B S -> (f 


>-) 


< 5.8 x 10" 8 


[43] 



Table 2: Flavor-changing observables used in the fit. The BR(£> — > X s £ + £ ) is intended in the 
range qj +i - G [1,6] GeV 2 . 



their current experimental determinations. Specifically, a quantitative test of the model is 
obtained through a \ 2 function defined as 



xm = — zt^ — . ( 3 -^ 



^ (^ 2 )exp + (^)thco ' 



where Oi indicates the experimental value of the observables and fi[d] the corresponding 
theoretical prediction, which will be function of the model parameters collectively indicated 
with i?. For each quantity, the experimental and theoretical standard deviations are also 
reported in tables 1 and 2. For those among the observables having a negligible exper- 
imental error, we took as overall uncertainly 0.5% of the experimental value, which we 
consider a realistic estimate of the numerical error associated with the calculations. For 
the details of the estimation of theoretical errors on the flavor observables (table 2), we 
refer the reader to the comments reported in ref. [19]. 

In evaluating the \ 2 function, we also included the bounds reported in tables 1 and 2, 
in the form of suitably smoothened step functions added to the x 2 function and returning 
zero in the case of a respected constraint. Concerning SUSY masses, the bounds used in 
our analysis are the most conservative ones reported by the PDG [26]. In particular, we do 
not use an explicit bound on the gluino mass, since that on the lightest chargino turns out 
to be strong enough. On the SUSY mass bounds used we will comment again in section 
4.3. 
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Sector 


# 


Parameters 


EauEe 

o o 


3 


an. Mr. 


SUSY-breaking 


6 


mi6, w-i/9, , wifj,, Aj7, ^4n 


SUSY (EW scale) 


2 


tan/3, /x 


neutrino 


1 




3rd generation Yukawa 


1 


yt = yb = y T = y» T 


light generation Yukawa 


6 


yu,c: yd, si ye,fi 


CKM 


4 


A, A, r) 



Table 3: Model parameters. Unless explicitly stated, they are intended at the GUT scale. 

The x 2 function is minimized using MIGRAD, which is part of the CERNlib library [44]. 
The minimization procedure guarantees the invariance of our conclusions under reparam- 
eterizations of the theory. 

The full set of free parameters describing our considered class of models is collected in 
table 3. For our purposes, we can take the grand unified group, on which we do not have 
to make any assumptions, to be broken to the SM group in one single step beneath the 
unification scale Mq. We allow for a percent level threshold correction 63 to the strong 
gauge coupling at the GUT scale. Our assumption of universal sfermion masses and split 
(but real) trilinear couplings amounts to six parameters in the soft SUSY breaking sector. 
The MSSM renormalization group equations (RGEs) for the GUT-scale parameters are 
solved between Mq and the EW scale, where we define [i and tan (3. 

In order to account, in this RG evolution, for the (possible) effects of right-handed 
neutrinos, present e.g. in SO(10) and required for the see-saw mechanism, we allow for 
the contribution of a third-generation neutrino Yukawa coupling (with initial condition 
yv T = yt) in all RGEs between Mq and the right-handed neutrino threshold Mr < Mq 
[45, 46, 47]. We stress that, since we are not considering any particular model construction 
and are only concerned with the phenomenological viability of Yukawa unification, it is 
sufficient in our setup to consider only a single right-handed neutrino threshold and the 
details of the neutrino Yukawa texture are not relevant. Therefore, we simply assume 
(Yv)ij = yt^i3^j3 at Mq. The concrete impact of the inclusion of Mr on our results will be 
discussed in section 4.4. 

It is worth stressing that the inclusion of the RH neutrino scale in the running as dis- 
cussed above allows to eliminate a potential source of large logarithmic GUT-scale threshold 
corrections to YU as well as to Higgs splitting. With the above said, it is difficult to exclude, 
in our fully general approach, the presence of additional GUT-scale threshold corrections 
to the Yukawa couplings. However, since our aim is to explore the predictive power of the 
hypothesis of YU, we will assume residual threshold corrections to be negligible, i.e. to 
leave YU a well defined hypothesis at the GUT scale. This assumption is realized in many 
concrete models, see e.g. [28]. 
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Our treatment of the RGE running, inclusion of threshold corrections, determination 
of the Higgs VEVs and calculation of the various observables is largely similar to that of 
Refs. [23, 19]. Here we just spell out a few improvements, namely: the use of two-loop 
RGEs for the soft sector, the Yukawa and the gauge couplings [48]; the inclusion of the 
full 3x3 flavor dependence of the Yukawa couplings (in place of neglecting the effect 
of Yukawa entries different from the 33 one). Concerning the latter point, while this in 
principle introduces 10 additional free parameters into the fitting procedure (cf. table 3), it 
does not pose a problem in practice (e.g. of fitting convergence). Indeed, the light Yukawa 
couplings and CKM entries are largely insensitive to changes in the SUSY parameters and, 
given also their hierarchical structure, can be determined very easily for each fit point. We 
further note that, while this approach is in principle valid for any Yukawa texture specified 
at Mq, we assumed exact YU to hold for the 33 elements of the Yukawa matrices in the 
basis where the down-type quark and charged lepton Yukawas are diagonal. In models 
predicting large 23 or 32 elements in Yukawa matrices, there would be small corrections to 
exact YU in this basis. 

Before concluding this section, we would like to add a further general comment on the 
procedure described above. It is clear that addressing the viability of the YU hypothesis 
is a non-trivial task also on purely computational grounds, since, in addition to the RG 
evolution of the relevant couplings, one has to take into account also threshold corrections 
to fermion masses, that depend on the details of the SUSY spectrum, as well as uncer- 
tainties in the measured values of fermion masses themselves and additional constraints 
on the parameter space, like FCNC processes. To deal with these problems, two different 
approaches have been followed in the literature. 

The first approach [22, 49, 50, 51, 17, 18, 52, 53, 54] amounts to fixing the low-scale 
values of the fermion masses to their experimental central values. After taking into account 
weak-scale threshold corrections and the running of the Yukawa couplings to the GUT 
scale, the differences between yt, y\, and y T as functions of the SUSY breaking parameters 
quantify the quality of the - approximate - YU. Apart from computational simplicity, this 
approach is motivated by the fact that exact YU might be spoiled by higher-dimensional 
operators and it allows to sacrifice some amount of unification to accommodate additional 
constraints, like the dark matter relic density. 

The second approach [27, 28, 23, 19], which we adhere to, amounts to imposing exact 
YU at the GUT scale, while fitting the low-energy values of fermion masses to their observed 
values by means of a x 2 minimization procedure. In our view, in spite of possibly being 
computationally more demanding, this approach has several advantages: First, the fitting 
procedure automatically singles out regions in parameter space favouring YU, while these 
regions have to be found by a widespread scan in the previous approach; second, the exact 
YU case represents a benchmark case exploiting the maximal predictive power of the YU 
hypothesis, hence leading to clear-cut predictions about qualitative features of the SUSY 
spectrum, which are univocally falsifiable at the LHC. This predictability has the downside 
that e.g. it may not be as easy to reconcile the standard neutralino relic density with YU 
in this scenario. However, we note that this constraint is a very indirect one, since it can 
easily be circumvented by departing from standard cosmology, without compromising the 
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Figure 1: Left: lines of constant \ 2 m the [i vs Ad plane, with mig = 4 TeV. Right: contributions 
to the x 2 function from BK(B — > X s 7). The black dot represents the example fit reported in 
tables 4-5. 



other ingredients of the approach. 
4. Numerical analysis 

In this section we would like to show that the trilinear splitting scenario, as defined at the 
end of section 2, allows generic SUSY GUTs with exact YU to recover full phenomeno- 
logical viability without invoking decoupling, at variance with what happens in the case 
of universal trilinears [23, 19]. This fact is due to a non-trivial interplay between the \i 
parameter and the Ajj — Ad splitting. The basic picture is that these parameters control 
the SUSY corrections to the bottom mass nib and to the B — > X s ^f decay rate, the latter 
representing the two main observables that generate a problem in the case of universal 
soft terms, once other constraints, notably BR(S S — > and BR(S — > X s £ + £~), are 

taken into account. In the numerical analysis, we restricted ourselves to a fixed universal 
sfermion mass miQ of 4 TeV, also for the sake of comparison with refs. [23, 19]. While 
a change in this parameter would change the allowed room for the stop-sbottom splitting 
mentioned in section 2, it would leave the lightest part of the spectrum, which is relevant 
in particular for collider phenomenology, unaffected as we will see. 

In the following subsection, we will present the main results of our numerical analysis, 
and their theoretical interpretation will be provided in section 4.2. Sections 4.3-4.5 are 
devoted to additional considerations, related to the implied sparticle mass spectrum, to 
the role of right-handed neutrinos, and to the (<? — 2) M and dark matter constraints. 

4.1 Main results 

The left-hand panel of fig. 1 displays the lines of constant x 2 ( as defined in eq. (3.1)) in 
the \x vs Ad plane. This plot is obtained by sampling the x 2 function on this plane via 
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MIGRAD minimizations where only mis, \i and are kept fixed. The right-hand panel of 
fig. 1 shows the contribution to the x 2 function solely from the BR(i? — > Xg'y) constraint, 
backdropped by the total x 2 contours for comparison. In fig. 2, we show the contours 
of the four most important input parameters of table 3 chosen by the fitting procedure. 
Finally, in figures 3-5, we show the contours of several masses of interest. 
Let us make some immediate observations. 

1. In the left-hand panel of fig. 1, the lowest values for the x 2 function are obtained for 
[i ~ 3.5 — 4.5 TeV, which is of the order of miQ, here set to 4 TeV. 

2. The x 2 contours are roughly symmetric with respect to the axis Ad = 0. 
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Figure 3: Lines of constant masses (in GeV) chosen by the fit for the lightest stop (left) and the 
gluino (right). 

3. From the top-left panel of fig. 2, one can see that the preferred region points to Ajj ~ 
—2.5 mi6, and in particular to a sizable A\j — Ad splitting. The bad phenomenological 
performance of YU in the case of universality between trilinears [23, 19] is recovered 
in the limit Ad — > Ajj « — 2m\Q, and /i <C mi6, i.e. in the lower, leftmost part of 
the plot. 

4. From the right-hand panel of fig. 1, one can see that the BR(i? — ► Xg'y) constraint is 
the main driving force for large fi. Note that, with our assumptions on the theoretical 
uncertainties as in table 2, xl^ sy = 1 (2) corresponds to BR(B -> X s7 ) ~ 3.0 x 10" 4 
(2.5 x 10~ 4 ). The branching ratio is always below the experimental central value due 
to the destructive interference between SM and chargino contributions. 

5. The top-right panel of fig. 2 shows that the universal gaugino mass m^j increases 
with increasing Ad, and is also dependent on the sign of Ad- For positive Ad and 
[i > 4 TeV, there is a steep rise in rn.i/2- 

6. For \i > 4.5 TeV, the x 2 starts deteriorating again, although the % 2 contribution 
from BR(S — > X s ^) is very small in this region. Both the stop mass and the gluino 
mass, as shown respectively in the left and right panels of fig. 3, increase again in 
this region. 

7. The two plots at the bottom of fig. 2 show that is preferred to be positive and 
equal to or less than mf 6 , while Tn/jj is preferred to be zero or even negative at the 
GUT scale. This is in contrast to refs. [23, 19], where both fn 2 H ^ and mjj were 
positive and greater than mf 6 . rn 2 Hu is basically independent of Ad- 

Points 1.-3. above confirm the previously obtained results [19] that YU is disfavoured 
for Ajj = Ad by an interplay between the corrections to the bottom quark mass and the 
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Table 4: Example of fit in the region with successful YU. The pull in the last column is defined 
as the square root of the x 2 contribution. 

FCNC constraints, but they also show that the trilinear splitting scenario considered here 
indeed gives rise to a viable solution featuring exact YU and being compatible with all 
relevant constraints. Interestingly enough, the recovery of phenomenological viability is 
not obtained by invoking a decoupling of the sparticle spectrum, but it instead seems to 
require parts of this spectrum to be very close to their experimental lower bounds. 

In table 4, we report the fitted values for the observables entering the y 2 function for 
one fit belonging to the region with lowest x 2 - The input values for this example fit are 
reported in the left panel of table 5 and the resulting spectrum predictions, on which we 
will comment again in section 4.3, on the right panel of the same table. This example fit 
is also represented in figures 1-3 as a dot. Apart from a pull in the Z mass, the largest 
contributions to the x 2 come from B — > X s "y and B — > rv. However, it should be noted 
that in both of these cases, the SUSY contribution at this parameter point is very small and 
most of the 'discrepancies' are rooted in discrepancies between the current experimental 
central values and the SM prediction. The prediction for B s — ► is well below the 

current experimental upper bound (cf. table 2) but well within the LHCb (and probably 
even Tevatron) reach. Since the predictions are in the (1 — 3) x 10 -8 ballpark in the 
entire preferred region, an experimental upper bound below 10 -8 would put this scenario 
in question. 

In the following section, we will attempt to interpret the above findings on theoretical 
grounds. 

4.2 Interpretation of fit results 

First of all, let us discuss the items 1.-3. of the above observations, i.e. the questions 
why the trilinear splitting helps to obtain successful YU, why this mechanism is roughly 
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symmetric under a sign change of Ad and why it requires large \i. The mechanism at work 
can be understood by recalling the basic ingredients for YU discussed in section 2, and it 
turns out to be quite compelling. 

Recall that, to suppress the 
gluino corrections to with re- 
spect to the chargino ones, a hier- 
archy mi. <C m~ h is required in ad- 
dition to a large trilinear parame- 
ter for the stop at the electroweak 
scale. This hierarchy is dependent 
on the trilinear parameters of the 
up- and down-type squarks, be- 
cause the latter contribute to the 
RG evolution of the running squark 
masses. Assuming tan (3 ~ 50 and 
neglecting gaugino mass contribu- 
tions (anticipating the condition 
mi/2 < "ii6, Au,d), the low-energy 
values for the third-generation squark 
mass-squared parameters can be 
approximately written in terms of 
the GUT-scale parameters as 
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Table 5: Input parameters and spectrum predictions for 
the example fit reported in table 4. All masses and massive 
input parameters arc in units of GcV. 



(m|) 33 « 0.51 m 2 l6 - 0.12 m 2 Hu - 0.09 m 2 Hd - 0.02 A 2 V - 0.02^ , 



( m D, 



0.49 mf 6 - 0.22 m 2 Hu - 0.01 m 2 Hd - 0.06 A% + 0.01 A D 

0.21 m 2 Hd +0.01 A% - 0.05 Zj, 



2 ^ ~ 0.55 m? 6 + 0.01 



(4.1) 
(4.2) 
(4.3) 



As is apparent from eq. (4.2), a very light right-handed stop can be obtained by appropri- 

ately adjusting m 2 H ^ and at the GUT scale. In the universal case, A v = A D = Aq, a 

sizable Aq also leads to a reduction of {rn 2 D ) 33 and (771^)33. Instead, in the non-universal 
— 2 — 2 

case, the choice A D <C A v allows to maintain a light right-handed stop, while preventing 
negative RGE contributions to the right-handed sbottom and left-handed squark masses. 
As a result, this mechanism permits to obtain a strong mass hierarchy 



(m 2 u) 33 « (™q) 33 < {m 2 D ) 



33 



(4.4) 



implying 7 



*fi *L b L b R > 



(4.5) 



which is what is needed to maximize the negative chargino corrections to nib and suppress 
the gluino contributions. While this hierarchy is also present in the universal case, it can 



7 Note that the LR mixing terms in the squark mass matrices do not play a role in this discussion, since 
mi 6 S> m t At and » mt,/xtan/3. 
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be greatly amplified in the trilinear splitting scenario by reducing A D : this leads to only 
a mild increase of and , but a strong increase of mg , while leaving th^ r almost 
unaffected. 

On the basis of equations (4.1)— (4.3), this RG effect is manifestly invariant under a sign 
change in either Ay or Aq. This invariance - which explains the approximate symmetry 
observed in item 2. - is in fact due to the condition <C A\j^, which is preferred by 
Yukawa unification, and would otherwise be spoiled by terms proportional to mi/2-Au,D in 
equations (4.1)-(4.3). 

We stress at this point that, in our numerical analysis, we calculate the physical mass of 
the light stop at the one-loop level, similarly to our procedure in ref. [19], since the lightness 
of the fitted stop mass implies that one-loop corrections are crucial to assess whether a given 
parameter point is viable or excluded by a tachyonic stop. For the remaining sparticles, 
we use the running masses. 

Once the above mechanism has ensured that the overall sign of the threshold correction 
to mj, is negative, an increase in \i, to which the threshold correction is proportional (see eqs. 
(2.2) and (2.3)) helps to make it parametrically large enough in magnitude, as required to 
fit the experimental data. At the same time, and quite interestingly, this large \i suppresses 
the chargino contributions to the b — > s-y amplitude (see [55] for a discussion on this point), 
therefore preventing a large destructive interference with the SM contribution. 8 This makes 
clear at the same time why the BR(i? — ► X s ^) constraint dominates the x 2 f° r small \x, as 
observed in point 4. of section 4.1 and shown in the right-hand panel of fig. 1. We stress 
that, similarly to ref. [19], we allow for both signs of the b — > S7 amplitude in our analysis. 
The large fi solution described above is preferred by the fit over the solution where the 
SUSY contributions are so large that they flip the sign of the b — ► 57 amplitude. In fact, 
the latter case turns out to imply a too large branching ratio of B — ► X s £ + £~ with respect 
to the experimental measurement [57]. 

Concerning point 5. of section 4.1, the value of m^ 2 chosen by the fit, this is due to the 
requirement of a light gluino mass to suppress the gluino corrections to mj,, as discussed 
in section 2. In fact, my 2 is always fitted close to its lowest allowed value, set by the 
LEP lower bound on the mass of the lightest chargino (see table 1), which is an almost 
pure Wino in our setup, due to the large //. If only one-loop RGEs for the gaugino masses 
were used, M2 > 104 GeV would imply m\/2 > 132 GeV. However, due to the conditions 
|^.£/,z>| ^ m i/2) two-loop effects become important in the running of the gaugino masses. 
These two-loop contributions are responsible both for the possibility of having m\/2 less 
than 132 GeV and for the rise of mjft with Ad, visible in the top- right panel of fig. 2. 
This effect always ensures a light chargino, except in the top-right corner of the plot. This 
region will be discussed in the next paragraph. 

For too large fx, the x 2 starts worsening again, as mentioned in point 6. of section 4.1: 
in fact, in this region, the negative corrections to nib start being so large that the mechanism 
above has to be tamed to prevent nib from dropping below 4.2 GeV. There are different 

8 For large fi x tan/3 one may worry about the size of gluino contributions as well [56, 55]. In our case, 
gluino contributions have roughly the same size as those from Higgses (in turn, an order of magnitude 
smaller than charginos), but opposite sign, thus cancelling with each other almost exactly. 
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possibilities to achieve this: for fi > 4.5 TeV and Ad < 0, the fits tune the lightest stop mass 
to be larger, as shown in the left panel of fig. 3, reducing the size of the chargino corrections 
to m&; for fi > 4.5 TeV and Ad > 0, the fits instead increase m 1( / 2 and accordingly the 
gluino mass, as shown in the right panel of fig. 3, increasing the gluino corrections to m^. 
While these mechanisms allow to obtain a correct value for nib, they cause tensions in other 
observables, leading to a steep rise in the x 2 ■ Therefore, viable YU solutions at even higher 
[i are not to be expected. 

Regarding the values of m 2 H commented on in point 7. of section 4.1, rrVjj^ is basically 
independent of Ad because its value is fixed by the EWSB conditions. Indeed, at large 
tan/?, these conditions require fn 2 H ^ ~ — |/i| 2 to hold at the EW scale. The value of m/jj 
on the other hand is bounded from below because EWSB requires rrvfj > m 'ji u a ^ t ne EW 
scale and bounded from above because a too large value would drive the sbottom masses 
smaller, cf. eq. (4.3), which is unfavourable for YU. 



The possibility to indeed fulfill the weak-scale conditions — 



m 



< 



m H , thereby 



achieving correct EWSB, can be illustrated through the following approximate expressions 
for the Higgs soft terms 



m 2 Hu {M z ) = —0.74 mf 6 + 0.56 rn? Hu + 0.06 m 2 Hd — 0.11 A v + 0.01 A D , (4.6) 
m 2 Hd {M z ) = -0.81 m? 6 + 0.06 m 2 Hu + 0.52 m 2 Hd + 0.01 A^j - 0.UA D , (4.7) 

where the parameters on the r.h.s. of either equation are at the GUT scale. These ex- 
pressions reproduce with remarkable accuracy the Higgs soft terms calculated with the full 
numerical procedure and may be used, e.g., on the parameter values of the example fit in 
table 5. 

For this fit, the evolution of the Higgs 
soft terms to low scales is also reported in 
figure 4. The latter shows in particular that 
fn 2 H d stays approximately constant through- 
out the running range, at variance with mjj , 
. io . / thus eventually allowing the weak-scale in- 

equality < m H d - Hence from this fig- 
ure and table 4 one can conclude that EWSB 
can indeed be quantitatively fulfilled through- 
out the explored parameter space, with pulls 
on EW observables approximately constant 
in the best-fit region. 9 

Considering that the natural EW scale 
is sensibly below the tuiq- and /i-ranges pre- 
ferred by the mechanism discussed above, it 
is clear that our scenario involves some de- 
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Figure 4: RG evolution of the parameters m 2 H 
for the example fit in tables 4-5. 



9 We note that inclusion of the one-loop tadpoles [58] in the EWSB equations plays in our case an 
important role, since these tadpoles will receive large logarithmic corrections from various among the heavy 
particles present in the spectrum. We thank the Referee for triggering this discussion. 
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Figure 5: Lines of constant masses (in GeV) chosen by the fit for the lightest sbottom (left) and 
the light neutral Higgs (right). 



gree of fine tuning in order for EWSB to be successful at the quantitative level. We also 
note that, in the case of the trilinear splitting scenario, the inverse scalar mass hierarchy, 
i.e. the hierarchy between the third and the first two generations sfermions, is reduced 
because of the heavier masses for the sbottoms and staus. This hierarchy, which improves 
the fine-tuning in the corrections to the Higgs mass, would instead be active in the case of 
universal trilinear couplings [27, 28, 23, 19]. We emphasize, on the other hand, that the 
trilinear splitting mechanism implies a very light stop mass, with a Higgs mass comfort- 
ably above the LEP bound, and with all FCNC constraints automatically fulfilled. On a 
model-dependent basis, all these desired features would require corresponding amounts of 
fine tuning as well, that in our case are simply absent. 

To summarize, the mechanisms described above single out a region in parameter space 
where successful YU is obtained, standing all the other experimental constraints, by a 
non-trivial interplay between the requirements of a large enough negative correction to mj 
and a small enough correction to b — > 57. As we saw above, allowing for split trilinear 
couplings, this can be achieved by \A^\ 2 <C \A\j\ 2 and large \x ~ m\% (but not too large). 
We emphasize that the recovery of phenomenological viability is not obtained by invoking 
a decoupling of the sparticle spectrum, it instead strongly requires parts of this spectrum 
to be very close to their experimental lower bounds. Since this observation is crucial for 
the LHC phenomenology of this class of models, we now briefly discuss the SUSY spectrum 
entailed by successful YU. 

4.3 Sparticle spectrum with successful YU 

As discussed in the previous sections, the region in parameter space favouring YU features 

10 We observe, however, that a certain amount of fine tuning in EWSB is by now a common feature within 
low-energy SUSY (for an insightful discussion, see [59]). 
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a gluino with a mass of around 350 GeV and a very light stop. In fact, as explained, a 
light stop mass is one of the main ingredients of the mechanism that allows successful YU 
in our trilinear splitting scenario. Therefore, a stop mass as low as about 100-200 GeV is 
generally preferred. However, as can bee seen in the left-hand panel of fig. 3, more generally 
stop masses < 400 GeV can lead to viable Yukawa unification. The large stop-sbottom 
splitting required by YU, as discussed in section 2, can be appreciated in fig. 5, showing the 
contours of the lightest sbottom mass, which is always nearly degenerate with the heavy 
stop, cf. (4.5). The lightest chargino and light neutralino masses are preferred to be as 
light as is allowed by experiment. In the right-hand panel of fig. 5, we show the contours 
of the light neutral Higgs boson mass. Interestingly enough, it is found to be in the 120- 
125 GeV region, such that the LEP constraint is not active. The heavy neutral, charged 
and pseudoscalar Higgs bosons are nearly degenerate and assume values around 1 TeV. All 
these predictions are exemplified in the right panel of table 5. Concerning the part of the 
SUSY spectrum not reported in this table, masses lie in the ballpark of mi§: for first and 
second generation sfermions, because of the absence of a large Yukawa contribution in the 
running; for the heavy neutralinos as well as the heavy chargino, because they are mostly 
Higgsinos, with masses of 0(/x). 

The non-decoupling and clear-cut nature of these spectrum predictions make the split 
trilinear scenario a compelling and testable one. In this respect, a first interesting question 
is whether part of our stop vs gluino parameter space may actually not comply with the 
SUSY mass bounds provided by Tevatron (see e.g. [60]), which are stronger than those 
in our table 1. We note however that the Tevatron bounds typically assume mSUGRA 
scenarios and we cannot identify any obvious relation to translate those bounds to our 
case. We believe that this issue would deserve a separate study. 

4.4 The role of right-handed neutrinos 

At this point, it is worth commenting on the effect of our inclusion of right-handed neutri- 
nos, as discussed in section 3, on the numerical results. 

Neutrino Yukawa couplings enter in the RGEs of up-type quark and charged lepton 
Yukawas and tend to drive these couplings to smaller values. In the leading-log approxima- 
tion, the difference between the values of the top and tau Yukawa couplings at low energies 
in the presence of neutrino Yukawa contributions and the values they would take in the 
absence of right-handed neutrinos is 



This percent level change in Yukawa couplings can be compensated by adjusting accordingly 
the GUT-scale value of the Yukawa coupling and tan /3, which in turn affects the value of 
the b quark Yukawa coupling at low energies. 

So, while right-handed neutrinos in principle affect the evolution of the third gener- 
ation Yukawa couplings, we found numerically that the presence of the neutrino Yukawa 
contributions and the value of the right-handed neutrino mass scale Mr do not signifi- 
cantly affect the success of YU for any given parameter point, since the small changes in 





- 19 - 



Yukawa couplings induced by these contributions can easily be compensated by changes in 
the remaining input parameters. 

In the soft sector, the neutrino Yukawa contributions leave the largest impact on the 
left-handed slepton doublet mass term m 2 L and on the up-type Higgs mass term rr? H . 
In the leading-log approximation and with the boundary conditions (2.11)-(2.12), their 
low-energy values are modified according to [45, 47] 

{Oij - K) -J = - ^ He + 2mk + 44) tfY^ log Qg) , (4.9) 

m 2 Hu - (m 2 Hu r = ( 4 ™?6 + 2m^ u + 4^) Tr{YjY v ) log (j^ , (4.10) 

where the quantities on the right-hand side are defined at the GUT scale. According to 

(4.9) , the presence of right-handed neutrinos leads to lighter left-handed sleptons at low 
energies; however, this does not have any relevant impact on the mechanism ensuring the 
success of YU. The off-diagonal components of (4.9) give rise to lepton flavor violating 
decays [45], but definite predictions can only be made in models predicting (yJYj,)^-. Eq. 

(4.10) shows that the neutrino Yukawa contributions drive mjj to smaller values. However, 
this can be easily compensated by raising the value of rn 2 ^^ at the GUT scale, which is 
possible in the setup of non-universal Higgs masses. We stress that this change in mjj 
induced by right-handed neutrino effects is not sufficient to explain or to generate the large 
m H u ~~ m # d splitting required for successful YU. 

To summarize, our approach of taking into account contributions from right-handed 
neutrinos on the evolution of couplings introduces one more free parameter, Mr, which 
allows to account for neutrino induced threshold corrections in particular to Yukawa cou- 
plings and to rn 2 H . However, numerically, the fits turn out to be quite insensitive to the 
value of this parameter and even a removal of right-handed neutrino effects by taking 
Mr — > Mq does not significantly affect the results. This means on the one hand that 
YU does not prefer or single out a particular scale for right-handed neutrinos; on the 
other hand, it means that the mechanism identified in our analysis cannot be spoiled by 
right-handed neutrino effects. 

4.5 Remarks on (g — 2) M and dark matter relic abundance 

Before concluding this section, we would like to comment on two observables we did not take 
into account in our numerical analysis: the muon anomalous magnetic moment a M = (<?— 2)^ 
and the neutralino cosmological relic density. 

Under our assumptions of sfermion mass universality, cf. (2.11), we found a posteriori 
predictions for the SUSY contributions to at the 10~ n level, much smaller than the 
current, O(10 -9 ) discrepancy between experiment and the SM prediction [61]. The reasons 
for this smallness are on the one hand the large values of the \i parameter and on the 
other the relatively heavy slepton spectrum, both leading to a suppression of the chargino 
contributions to a«. We note however that a relaxation of the universality assumption 

(2.11) and the inclusion of a M in the fitting procedure could help reach the 10 -9 level, 
while hardly affecting the mechanism discussed in section 4.2. We observe nevertheless 
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that the theoretical and experimental status of the {g — 2)^ tension still remains to be 
settled. 

Concerning the WMAP constraint on the dark matter relic density, again we did not 
include this observable in our numerical analysis, since we consider it to be a very indirect 
constraint which can be spoiled by many cosmological effects. A posteriori, and under the 
assumption of a completely "standard" thermal history of the universe, we found values 
of the relic density much higher than allowed. However, it may be possible to sufficiently 
suppress the predicted relic density in the case it were included in the fitting procedure. 
In fact, since the lightest neutralino is always at the level of roughly 60 GeV, it might be 
possible, albeit through some fine tuning, to exploit the Higgs funnel region occurring at 
2m x o ~ rah. For a detailed study of this possibility, see ref. [52]. 

5. Model-building discussion 

The pattern of soft SUSY-breaking terms considered in this paper points to mechanisms 
of SUSY-breaking that do not, in general, respect the GUT group, and that, on the other 
hand, have a highly specific flavor structure. For example, in the general case of GUT-scale 
MFV, the only flavor spurions generated should be the SM Yukawas. These mechanisms 
would be at work at or above the unification scale Mq- 

As discussed in section 2, the trilinear splitting scenario of eqs. (2.11)-(2.12) analysed 
in this paper is a special case of MFV soft terms. As illustrated in fig. 2, for miQ = 4 TeV, 
our fitting procedure to low-energy data points to the following patterns of soft-breaking 
terms: 

Au -2.5 mi 6 , < \A D \ < \A V \ , 

m 2 Hu , \m 2 Hd \ = 0(m? 6 ) , (5.1) 
m 1/2 <. m 16 . 

In this section we would like to discuss a concrete example of a SUSY-breaking scenario 
where this pattern is naturally realized. Our example aims at illustrating to which extent, 
in a very predictive framework like that of Yukawa unification, the combined information 
from existing low-energy data can be translated into information on the mechanism of 
SUSY breaking at work. 

The one pattern most clearly emerging from our analysis is a simultaneous splitting 
in nig — rn? Hd and in Au — Ar>. This hints at SUSY-breaking VEVs coupled to operators 
that distinguish the 'up' and 'down' directions in Higgs soft-terms and ^4-terms alike. The 
arguably simplest scenario where this can be realized within SO(10) is with a spurion field 
that gets a VEV in the adjoint representation of SU(2)#, namely V ~ (1, 1,3) (recall, in 
this respect, that SO(10) contains the Pati-Salam group SU(4) x SU(2) L x SU(2) B [62]). 
One can further assume the presence of a spurion field X, that is instead a singlet under 
SO(10). A linear combination of the X and V fields will then introduce F-term SUSY 
breaking of the form c\Fx + C2-Fy, with Fy proportional to the T3 generator of SU(2)#, 
distinguishing the 'top' from the 'bottom' direction. 
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A-term splitting will then be simply generated by the superpotential term Wa oc 
Q{c\X + C2V)HQ, where Q is the left-handed Pati-Salam chiral supermultiplet, trans- 
forming as Q ~ (4,2,1), and containing the MSSM chiral supermultiplets Q and L (see 
[63] for an introduction to the formalism). The field H contains instead the MSSM H u 
and Hd superfields and transforms as H ~ (1,2,2). After the breakdown of SUSY, the 
coupling Wa will then generate ^4-terms of the form Ayp oc c\(Fx) ±C2(Fy) /2, the overall 
coefficient being fixed by the Wa normalization. 

The above field content easily gives rise to split m 2 H ^-terms as well. In fact, operators 

of the form Kjj oc X1X2H* H , will naturally be present in the Kahler potential, where the 
fields Xi can be chosen as X or V, with gauge indexes adjusted so as to satisfy invariance 
under the Pati-Salam group. Among these operators, those inducing a single power of Fy 
will distinguish H u from Hj, thereby splitting the corresponding soft terms. 

It is worth observing that all the mass scales generated through the described mecha- 
nism are of the order |(^v;x) 2 |/Mp laIlck) which, in turn, will be of the order the gravitino 
mass, namely the TeV scale. 

The above mechanism will also, in general, produce soft terms for squark and slepton 
bilinears. In particular, there is no obvious symmetry argument by which the Kahler po- 
tential operators mentioned above for Higgs bilinears would not induce sfermion bilinear 
splittings as well. In our paper, sfermion bilinears have been assumed all degenerate to the 
value toiq. The absence of splittings in our case has been justified on purely phenomeno- 
logical grounds, namely, within the trilinear splitting scenario, data do not require bilinear 
splittings as well. The possibility that the above splitting mechanism could practically be 
ineffective for squark and slepton bilinears may be justified in frameworks where SUSY is 
broken through orbifold compactifications, and Higgs superfields live in the bulk whereas 
quark and lepton superfields are localized on the PS brane and the spurions X and V live 
on the SO(10) brane (see [64] for an example). Finally, since mi/2 is generated from the 
gauge kinetic function, it is completely unrelated to the above mechanism and can well be 
smaller in magnitude than mig, as in eq. (5.1). 

The discussion in this section also highlights the importance of information on the 
lightest part of the SUSY spectrum, e.g. from the LHC. If the pattern were compatible 
with that described in section 4.3, this would allow, depending on, say, M g vs (see 
respectively figs. 3, right and 5, left) to virtually select a point in the A& vs fi plane (see 
the panels of fig. 2). The implied information on the main input parameters would permit 
correspondingly sharper model-building considerations than those presented above. 

6. Conclusion 

We have considered general SUSY GUT frameworks with exact Yukawa unification and 
where the hypothesis of universal GUT-scale soft terms is relaxed. We have first entertained 
the general possibility that soft terms be of minimally flavor violating form. In this case 
the hypothesis of exact YU and the hierarchical structure of the Yukawa couplings allow 
to parameterize squark bilinears and trilinears in a general, simple and accurate way as in 
eqs. (2.9)-(2.10). Among the soft-term non-universalities allowed by this general param- 
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eterization, we have then focused on the scenario where up-type and down-type trilinear 
soft terms are split from each other. 

We have explored the viability of this trilinear splitting scenario by contrasting the 
model predictions for EW observables, quark masses and quark FCNC processes against 
data in a global fitting procedure. Agreement with data singles out one main scenario, 
featuring a sizable splitting between the j4-terms and a large //-term. In spite of a slight 
increase in the fine tuning required to achieve EWSB with precisely the correct value 
of Mz, this scenario allows a substantial improvement on other observables that, on a 
model-dependent basis, do often require some amount of fine tuning as well. First, and 
quite remarkably, phenomenological viability does not invoke a partial decoupling of the 
sparticle spectrum, as in the case of universal soft terms, but instead it requires part of the 
spectrum, notably the lightest stop, the gluino and the lightest chargino and neutralinos, 
to be very close to the current experimental limits. The lightest Higgs particle is also well 
above the LEP bound, it is actually quite robustly predicted at around 125 GeV. Second, 
the above parameter space is selected by a non-trivial interplay between the requirement of 
negative, sizable SUSY threshold corrections to nib, and an instead negligible modification 
of the B — > X s "f decay rate, in presence of various other EW and -B-physics constraints. 
Hence the very same mechanism that makes the correction large enough, automatically 
allows all FCNC constraints to be fulfilled. 

We have also discussed a possible model of SUSY breaking where the pattern of soft 
terms, selected above on sheer phenomenological grounds, is realized. This discussion 
highlights the crucial role of SUSY spectrum determinations at the LHC for either falsifying 
YU or else offering important hints on the mechanism of SUSY breaking at work. 

Our results provide a concrete example where, exploiting the full predictive power of 
YU and under what we consider very plausible assumptions for soft terms, enough remnant 
information on the high-energy symmetries does indeed survive at low energies for these 
symmetries to be reconstructible. They open up two natural directions of investigation. 
On the one side, a more complete exploration of the above discussed general pattern of 
GUT-scale soft term non-universalities and its connection with a plausible mechanism of 
SUSY breaking. On the other side, a closer look at the collider signatures this pattern 
points to. Both directions are subject of future work. 
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